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ABSTRACT 


An  analytical  and  experimental  investigation  is  made  of 
the  non-linear,  large  amplitude,  high  angle  of  attack,  stall 
flutter  behavior  of  cantilevered  graphite/epoxy  wings.  Ten 
six-ply  graphite/epoxy  wings  with  a  wide  range  of  bending- 
torsion  characteristics  were  constructed  and  styrofoam  fair¬ 
ings  epoxied  to  these  to  form  NACA-0012  airfoil  shapes.  Wind 
tunnel  tests  on  these  cantilevered  wings  revealed  torsional 
and  bending  stall  flutter  limit  cycles,  depending  on  the 
layup.  Reasonable  agreement  with  steady,  non-linear  theory 
and  with  unsteady,  linear  theory  was  found.  Fourier  analysis 
applied  to  the  ONERA  2-dimensional,  non-linear,  unsteady, 
aerodynamic  model  shows  reasonable  agreement  with 
2-dimensional  experiments  on  aerodynamic  force  and  moment 
coefficients.  A  3-dimensional,  non-linear,  unsteady,  aero- 
elastic  analysis  with  a  Newton-Raphson  solver  applied  to  the 
harmonic  balance  method  is  developed  to  attempt  to  predict 
the  non-linear,  stall  flutter  observed  in  the  experiment. 
Final  results  of  this  analysis  are  currently  being  completed, 
and  will  appear  shortly  in  a  forthcoming  report. 
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FOREWORD 


This  report  describes  work  done  at  the  Technology  Labo¬ 
ratory  for  Advanced  Composites  (TELAC)  at  the  Massachusetts 
Institute  of  Technology  for  the  Air  Force  under  Contract 
No.  F49620-86-C-0066 .  Dr.  Anthony  K.  Amos  was  the  contract 
monitor. 

The  work  reported  herein  was  performed  during  the  period 
1  July  1987  through  30  June  1988.  The  work  represents  the 
efforts  of  one  graduate  student,  Peter  Dunn,  and  an  under¬ 
graduate  student,  Chakko  Kovoor,  under  the  direction  of  Prin¬ 
cipal  Investigator,  John  Dugundji,  and  the  supporting  labora¬ 
tory  staff. 
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1 ■  INTRODUCTION 


The  present  research  is  part  of  a  continuing  investiga¬ 
tion  into  the  aeroelastic  flutter  and  divergence  behavior  of 
forward-swept,  graphite/epoxy  composite  wing  aircraft.  The 
specific  objectives  of  the  current  investigation  are  to 
investigate  experimentally  and  analytically,  the  roles  of 
non-linear  structures  and  non-linear  aerodynamics  in  large 
amplitude,  high  angle-of-attack  stall  flutter  of  aeroelasti- 
cally  tailored  wings. 

In  previous  investigations  at  M.I.T.,  the  aeroelastic 
flutter  and  divergence  behavior  of  cantilevered,  unswept  and 
swept,  graphite/epoxy  wings  was  investigated  in  a  small,  low- 
speed  wind  tunnel.  The  wings  were  six-ply,  graphite/epoxy 
plates  with  strong  bending-twisting  coupling.  Experiments 
were  conducted  to  determine  the  flutter  boundaries  of  these 
wings  both  at  low  and  high  angles  of  attack,  stall  flutter 
often  being  observed  in  the  latter.  Steady,  non-linear  aero¬ 
dynamics  correlated  well  before  the  onset  of  flutter,  and  the 
divergence  and  flutter  results  at  low  angles  of  attack  corre¬ 
lated  well  with  linear,  unsteady  theory,  indicating  some  ben¬ 
eficial  effects  of  ply  orientation  in  aeroelastic  behavior 
(References  1,  2,  and  3) .  But  no  non-linear,  unsteady  analy¬ 
ses  were  attempted  for  the  higher  angles  of  attack. 

Recently,  Tran  and  Petot  (Reference  4)  and  Dat  and  Tran 
(Reference  5)  of  Office  National  d’Etudes  et  de  Recherches 
Aerospatiales  have  developed  a  semi-empirical,  unsteady,  non¬ 
linear  model  (called  the  ONERA  model)  for  determining  the 
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2-dimensional  aerodynamic  forces  on  an  airfoil  oscillating  in 
pitch  only  which  experiences  dynamic  stall.  This  model 
incorporates  a  single  lag  term  operating  on  the  linear  part 
of  the  airfoil's  static  lift  curve,  thus  analogous  to  the 
Theodorsen  function  for  linear,  flat-plate  theory,  and  a  two 
lag  term  operating  on  the  non-linear  (i.e.  stalling)  part  of 
the  airfoil’s  static  lift  curve.  The  semi-empirical  coeffi¬ 
cients  of  the  non-linear  aerodynamics  for  the  ONERA  model 
were  determined  for  various  airfoils  and  the  model  applied  as 
a  comparison  against  experiment  by  Dat,  Tran,  and  Petot. 
Further  analysis  of  the  model  was  done  by  Peters 
(Reference  6)  to  differentiate  the  roles  of  angle  of  attack 
due  to  pitching  (0)  and  angle  of  attack  due  to  plunging 
(h/U) ,  and  by  Petot  and  Dat  (Reference  7)  to  reformulate  the 
equations  so  that  they  reduce  to  the  Theodorsen  function  in 
the  case  of  a  flat  plate  in  the  linear  domain.  Additionally, 
Petot  and  Loiseau  (Reference  8)  have  contributed  corrections 
to  the  ONERA  model  for  low  Reynold’s  number  flows,  in  the 
regime  of  the  current  investigation.  Generally,  however, 
little  work  has  been  done  in  implementing  the  ONERA  model  in 
a  non-linear,  aeroelastic  flutter  analysis.  The  present 
research  plans  to  use  the  non-linear  ONERA  model  to  help  ana¬ 
lytically  explain  some  of  the  stall  flutter  behavior  observed 
earlier  in  the  experimental  tests. 


5 


The  present  investigation  deals  with  the  effects  of 
large  amplitude  deflections  and  non-linear,  unsteady,  stalled 
aerodynamics  on  the  aeroelastic  behavior  of  aeroelastically 
tailored  wings.  Six-ply,  graphite/epoxy  wings  with  a  wide 
range  of  bending-torsion  characteristics  were  previously  con¬ 
structed  with  styrofoam  fairings  epoxied  to  these  to  form 
NACA-0012  airfoil  shapes.  These  wings  were  mounted  verti¬ 
cally  with  cantilevered  root  conditions  inside  the  M.I.T. 
Department  of  Aeronautics  and  Astronautics  5x7  ft  (1.5  mx 
3  m)  low  speed  acoustic  wind  tunnel,  and  tested  to  velocities 
of  30  m/s.  Figures  1  and  2  show  the  wing  construction  and 
the  typical  setup  in  the  wind  tunnel.  Ten  different  ply 
layups  were  tested,  namely  [02/90]g,  [90/02]g,  [+152/0]s, 
[-152/0]s,  [+302/0]s,  [-302/0]s,  [±15/0 ] s ,  [T15/0]s, 
[±30/0] _,  and  [*30/0]  .  These  wing  surfaces  include  some  of 

S  w 

those  used  in  previous  cantilever  tests  (References  1,  2, 
and  3) ,  so  that  the  reproducibility  of  those  previous  tests 
could  be  verified. 

Static  bending  and  free  vibration  bench  tests  were  con¬ 
ducted  to  verify  the  mass,  stiffness,  and  bending-torsion 
properties  of  the  wings.  The  wind  tunnel  tests  included 
divergence  and  limit  cycle  flutter  testing  at  both  low  and 
high  speeds,  with  root  angle  of  attack  varying  from  zero  to 
more  than  twice  the  static  stall  angle.  Data  was  recorded 
for  wing  bending  and  wing  torsion  moments  on  floppy  disk  via 
a  Nicolet  digital  oscilloscope  using  strain  gauges  located 


near  the  roots  of  the  wings.  Video  movies  were  also  taken  to 
record  tip  deflection  and  tip  twist.  The  experimental  data 
from  these  test  was  processed  and  analyzed  during  the  current 
reporting  period,  A  typical  example  of  some  of  the  observed 
flutter  boundaries  is  indicated  in  Figure  3. 

An  analytical  investigation  was  also  made  concurrently 
with  the  experimental  tests.  This  involved  a  5-degree  of 
freedom  Rayleigh-Ritz  flutter  analysis  using  the  5  modes, 
qji  (t)  ,  similar  to  those  of  References  2  and  3  (i.e.  1st  and 
2nd  wing  bending,  1st  and  2nd  wing  torsion,  and  1st  chordwise 
bending).  The  wing  torsion  modes,  which  previously  only 
reflected  the  root  warping  effects  on  the  natural  frequencies 
of  the  wings,  were  altered  to  also  reflect  the  root  warping 
effects  on  the  tip  deflection,  and  were  also  modified  to 
include  a  simple,  non-linear,  cubic  stiffening  effect. 

The  aerodynamic  analysis  employed  a  2-dimensional, 
unsteady,  non-linear  strip  theory  with  spanwise  correction 
for  finite  span  effects.  The  aerodynamic  forces,  namely  the 
lift  coefficient  CL,  the  moment  coefficient  CM,  and  the  drag 
coefficient  CD,  were  taken  from  the  ONERA  model  (Reference  7) 
in  the  form: 


Czi  *  szct  +  kV20  + 

CYl  +  X2CYl  -  Xz  [aozcc  +  oz0  ]  +  az  [aoza  +  oz*0  ] 
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C z2  +  2dwCz2  +  W  *1+ci  *  Cz2 


=  -w2(l+d2) 


9AC. 


+ei>rU 


where. 


a  =  0  - 

*  a 

(  )  2  - 


* 

h 


{  ) 

at 


and. 


T 


Ut 
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0  =  instantaneous  angle  of  attack 

h  =  instantaneous  deflection  of  1/4-chord 

h  *  =  non-dimensional  deflection 

b 

a  =  effective  angle  of  attack 

and  where  aQZ  is  the  slope  of  the  linear  part  of  the  static 
force  curve,  ACZ  is  the  non-linear  deviation  from  the 
extended  linear  lift  curve,  and  s2,  kV2,  Oz,  Xz,  5Z,  d,  w, 
and  e  are  previously  determined,  semi-empirical  lag  coeffi¬ 
cients.  In  the  above  formulation,  Czl  represents  the  usual 
linear  force,  while  Cz2  represents  the  additional  non-linear 
stalling  contribution.  The  Cy.^  term  is  the  circulatory  part 
of  the  linear  lift  force.  Fourier  analysis  is  used  to 
extract  the  low  harmonic  components  (1st  and  2nd  harmonics) 
from  the  non-linear  forcing  terms  of  these  differential  equa¬ 
tions,  by  assuming  simple  harmonic  motion  for  h(t)  and  0(t) . 

A  harmonic  balance  method  is  then  used  to  set  up  the 
coupled,  non-linear  equations  of  motion,  which  incorporate 
the  force  coefficients  from  the  above  equations.  First,  the 
components  of  the  force  coefficients  are  determined  assuming 
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a  sinusoidally  varying  input  for  the  angle  of  attack  and  the 

1/ 4-chord  deflection, 


0(X)  =  0O  +  0ssin(kx)  +  0ccos  (kx) 
h  (X)  *  h0  +  hssin(kX)  +  hrcos (kx) 

The  circulatory  part  of  the  linear  contribution  to  the 
2-dimensional  aerodynamics  is  given  by, 

Oy10U)  “  aoz®o 

CYls(x)  "  F(k)La(x)  -  G(k)Lc(x) 


Cylc(x)  «  G(k)Ls(x)  +  F(k)Lc(x) 

where,  in  the  present  analysis,  the  F  and  G  functions  are  the 
resulting  single  lag  approximations  to  the  Theodorsen  func¬ 
tion,  C (k)  «  F(k)+iG(k),  namely. 


F  (k) 


.2  .2 
kz  +  a2k 

xl  +  k2 


G  (k) 


Xzk<l-az) 
^  +  k2 


and  where  the  intermediate  variables,  Lg  (x)  and  Lc (x) ,  are 
determined  from  0Q,  0g,  8C  and  hD,  hg,  hc. 


La(x)  -  aOz[0s  +  khc  (x)]  -  C2k0c (x) 
Lc(x)  -  aO2[0c(x)  -  khs  (x) ]  +  Czk8s(x) 
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is  then 


The  complete  linear  contribution  C  lo,  Czlg,  Czlc 

found  by  adding  in  the  average,  sine,  and  cosine  components 

★  ★  * 

of  the  non-circulatory  terms  s2<X  and  kV2  &  . 

Similarly,  the  components  of  the  non-linear  contri¬ 
bution,  cz2o'  Cz2s'  anc*  Cz2c'  are  determined  from  the  ACzo, 
ACZS,  and  ACzc  components  of  the  non-linear  deviation, 
AC2  (T)  ,  as, 


The  equations  of  motion  for  the  bending-torsion  flutter 
of  the  wings  were  set  up  in  harmonic  balance  form  using  the 
5  modes,  q^ (T)  ,  mentioned  earlier,  and  considering  only  the 
constant  part  and  lowest  harmonic  of  each  coordinate, 

q.  (X)  =  q.  +  q.  sin(kx)  +  q,  cos  (kX) 

1  10  IS  1C 
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The  resulting  form  of  the  equations  of  motion  is  given  by  the 
3n  coupled  equations: 


"  [K] 

0 

0 

"{Q0r 

0 

-G)2  [M]  +  [K] 

0 

{qs} 

= 

{Qs} 

0 

0 

-to2  [M]  +  [K]_ 

_{qc>_ 

where  {qQ},  (qs),  anci  are  vectors  representing  the  har¬ 
monic  components  of  the  five  coordinates,  q^  (X)  ,  and  [M]  and 
[K]  are  the  appropriate  mass  and  stiffness  terms.  The  { Qq } , 
{Qs },  and  {Qc)  are  the  corresponding  harmonic  modal  forces 
obtained  by  integrating  the  non-linear,  2-dimensional  air 
forces,  Cz,  over  the  span  of  the  wing,  after  having  expressed 
the  angle  of  attack  0  and  deflection  h  in  terms  of  the 
coordinates  q^ (X)  . 

The  non-linear  equations  of  motion  are  then  solved  using 
a  Newton-Raphson  scheme  similar  to  that  used  by  Kuo,  Morino, 
and  Dugundji  (Reference  9)  for  a  related,  non-linear  panel 
flutter  problem.  The  3n  non-linear  equations  are  expressed  in 
the  form  f(x)=0,  and  the  sine  component  of  one  mode  q^ (x)  is 
set  to  some  small  constant  to  set  the  amplitude  level,  while 
its  cosine  component  is  set  to  zero,  since  the  flutter  limit 
cycle  oscillations  can  start  at  an  arbitrary  phase.  The 
Newton-Raphson  scheme  then  uses  the  3n  equations  of  motion  to 
solve  for  the  remaining  3n-2  modal  component  amplitudes,  and 
for  the  frequency  of  oscillation,  k,  and  flutter  velocity,  V. 
Because  the  sinusoidal  component  of  one  mode  shape  has 
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already  been  set  to  a  non-zero  value,  the  Newton-Raphson 
scheme  does  not  converge  to  the  trivial  steady  solution.  The 
Newton-Raphson  scheme  takes  an  initial  guess  of  the  state 
vector  x  and  drives  the  residual  of  the  equations  f (x)  toward 
zero  by  inverting  the  Jacobian  matrix  (derivative  matrix)  and 
obtaining  a  correction  Ax  to  the  current  guess.  The  process 
is  repeated  until  the  correction  Ax  becomes  negligible. 


Ax 


(n) 


df 


dx 


-l 


f  <x(n)) 


Jn 


(n+1) 


(n) 


+  Ax 


(n) 


This  Newton-Raphson,  harmonic  balance  scheme  has  already  been 
used  to  obtain  the  linear  flutter  solution,  and  is  currently 
being  applied  to  obtain  the  corresponding  non-linear  flutter 
solution  by  working  up  from  the  linear  solution. 

Comparison  of  experimental  flutter  boundaries  with  those 
obtained  from  previous  experiments  indicates  good  agreement, 
thus  demonstrating  good  reproducibility  of  the  experiments. 
The  static  tip  deflection  and  static  tip  twist  for  all  the 
wings  were  reasonably  predicted  by  the  steady,  non-linear 
analysis,  as  indicated  in  Figures  4  and  5.  The  unsteady, 
linear  analysis  reasonably  predicted  the  experimental  flutter 
boundaries  at  low  angles  of  attack.  The  Fourier  decomposi¬ 
tion  used  here  for  the  ONERA  model  reasonably  predicted  the 
2-dimensional,  unsteady  experiments  of  McAlister,  Pucci, 
McCroskey,  and  Carr  (Reference  10),  as  indicated  in  Figures  6 
and  7.  The  unsteady,  non-linear  analysis  described  here  is 
almost  completed.  The  trends  in  the  experimental  data  indi- 
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cate  that  the  average  unsteady  tip  deflection  and  tip  twist 
depart  from  the  predicted  static  tip  deflection  and  static 
tip  twist  once  flutter  oscillations  begin,  as  indicated  in 
Figures  4  and  5.  These  results  appear  to  indicate  a  non¬ 
linear  coupling  between  the  non-linear,  static  divergence 
problem  and  the  non-linear,  stall  flutter  problem. 

Final  results  of  the  present  investigation  are  currently 
being  completed  and  will  appear  in  a  M.S.  thesis  by  the  first 
author,  in  the  near  future.  This  forthcoming  thesis  will 
present  the  final  results  of  the  unsteady,  non-linear  analy¬ 
sis  and  will  compare  those  results  in  detail  with  the  experi¬ 
mental  data. 

IL  .ACCOMPLISHMENTS 

An  experimental  and  analytical  investigation  is  being 
made  of  the  linear  and  non-linear  aeroelastic  flutter  and 
divergence  behavior  of  graphite/epoxy  cantilevered  wings. 
Experimental  wind  tunnel  flutter  data  was  obtained  at  low  and 
high  angles  of  attack  in  a  low  speed  wind  tunnel,  for  can¬ 
tilevered  wings  with  NACA-0012  airfoil  shapes,  and  with  a 
wide  range  of  bending-torsion  characteristics. 

Stall  flutter  limit  cycles  were  observed  experimentally 
with  strain  gauge  and  video  data  taken.  Torsional  and  bend¬ 
ing  flutter  were  observed  depending  on  layup.  These  experi¬ 
ments  confirm  earlier  linear  analytical  findings  and  indicate 
a  non-linear  coupling  of  oscillatory  and  non-oscillatory 
aerodynamics  resulting  in  deviation  from  the  predicted  steady 
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analysis.  An  analytic  method  has  been  developed  to  incorpo¬ 
rate  non-linear  aerodynamics  into  the  aeroelastic  model,  and 
to  solve  the  resulting  coupled,  non-linear  flutter  problem. 
In  the  current  case,  the  non-linear  flutter  problem  couples 
with  the  non-linear,  stalled  divergence  problem,  and  the  two 
problems  must  be  solved  together. 

The  present  investigation  enlarges  the  experimental  base 
for  aeroelastic  tailoring  with  composites,  extending  the  data 
base  into  the  range  of  non-linear,  high  angle-of-attack 
flutter.  Along  with  the  corresponding  theoretical  analyses, 
these  should  provide  insight  into  the  actual  aeroelastic 
behavior  of  composite  wings  experiencing  dynamic  stall . 
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Figure  1  -  Wing  construction  and  dimensions 
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FLUTTER  VELOCITY  (M/S) 
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STALL  FLUTTER  EXPERIMENT 

FLUTTER  BOUNDARIES  OF  NACA-0012  WINGS 


Figure  3  -  Experimental  flutter  boundaries  for  [02/90]g, 
[+152/0]s,  [-152/0Js,  and  [-302/0]g  NACA-0012  wings 
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[90/09]  NACA-0012  WING,  NO  SWEEP 


[90/0-]  NACA-0012  WING,  NO  SWEEP 


(93G)  XSIMi  a 


ROOT  ANGLE  OF  ATTACK  (DEG) 


Figure  6  -  Hysteresis  loops  for  2-dimensional  lift 
coefficient  versus  angle  of  attack/  NACA-0012  air¬ 
foil,  Re  »  0.49  E6  (corrected  for  wind  tunnel  ef¬ 
fects)  ,  compared  against  McAlister,  Pucci,  McCroskey, 
and  Carr  (Reference  10)  .  - Analysis  - Experiment 
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THETA  (dag)  THETA  (dag) 


9?  dtgt 


Figure  7  -  Hysteresis  loops  for  2-dimensional  moment 
coefficient  versus  angle  of  attack,  NACA-0012  air¬ 
foil,  Re  *  0.49  E6  (corrected  for  wind  tunnel  ef¬ 
fects)  ,  compared  against  McAlister,  Pucci,  McCroskey, 
and  Carr  (Reference  10) .  - Analysis  - Experiment 
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